%simyo

p3_min = 3;
p5_min = 5;

p3_mega = 3;
p5_mega = 3;

megas3 = 500;
megas5 = 300;

cuota3 = 700;
cuota5 = 0;

c3_min = 0;
c5_min = 700;

minutos = 1:1:1000;
megas = 1:1:1000;

g5 = zeros(length(minutos), length(megas));
g3 = zeros(length(minutos), length(megas));
%Gasto por minuto
for i=1:length(minutos)
    for j = 1:length(megas)
        g5(i,j) = cuota5 + p5_min*minutos(i);
        g3(i,j) = cuota3 + p3_min*minutos(i);
        if (megas(j)>300)
            g5(i,j) = g5(i,j) + (megas(j)-300)*p5_mega;
        end
        if(megas(j)>500)
            g3(i,j) = g3(i,j) + (megas(j)-500)*p3_mega;
        end
        if (g5(i,j)<c5_min)
            g5(i,j) = c5_min;
        end
        if (g3(i,j)<c3_min)
            g3(i,j) = c3_min;
        end
    end
    
end
% 
% 
% subplot(131)
% title('<300mb')
% xlabel('minutos')
% ylabel('centimos')
% hold on
% plot(minutos, g3_300mb,'r');
% plot(minutos, g5_300mb,'b');
% hold off
% subplot(132)
% title('300mb-500mb')
% plot(minutos, c3 + t3_m*minutos,'r');
figure(1)
subplot(131)
mesh(minutos, megas, g3);
subplot(132)
mesh(minutos, megas, g5);
subplot(133)
mesh(minutos, megas, g3-g5);